Joint description of the e + e annihilation into both four-pion channels 
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The e + e~ — > 7r + 7r~7r°7r° reaction cross section as a function of the incident energy is calculated 
using a model that is an extension of our recently published model of the e + e~ annihilation into four 
charged pions. The latter considered the intermediate states with the n, p, and ai mesons and fixed 
the mixing angle of the a\p-K Lagrangian and other parameters by fitting the cross section data. Here 
we supplement the original intermediate states with those containing oj(782) and h\(\ f70), but keep 
unchanged the values of those parameters that enter both charged and mixed channel calculations. 
The inclusion of uj is vital for obtaining a good fit to the cross section data, while the intermediate 
states with hi further improve it. Finally, we merge our models of the e + e _ — > 7r + 7r~7r°7r° and 
e + e~ — > 7r + 7r~7r + 7r - reactions and obtain a simultaneous good fit. 
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The electron-positron annihilation into four pions has 
been theoretically studied by several authors [l|, [H, H, 0, 
5]. Assuming the one-photon approximation and vector 
meson dominance, this process goes via the p(770) meson 
and its recurrences. If some conditions on the rho-decay 
amplitude into four pions are met [(J , the cross section of 
the e + e~ annihilation into four pions at invariant energy 
W can be expressed in terms of the decay width of a p 
meson with mass W into four pions [!, 0|- Some of the 
models of the four-pion decay of the p meson H H, |J 
[To| can thus be conveniently utilized when determining 
the excitation function of the e + e~ annihilation into four 
pions. 

In our recent work Q we calculated the excitation 
function of the e + e~ annihilation into four charged pions 
and compared it to the existing data. We confirmed the 
conclusion of several experimental and theoretical papers 
@; EL IE HH Q3 that the axial- vector isovector resonance 
ai(1260) plays an important role. The new feature of our 
approach was that we did not take some a priori cho- 
sen Lagrangian of the aipir interaction, but considered 
a two-component Lagrangian that contained two param- 
eters: a mixing angle and an overall coupling constant. 
We varied the mixing angle in an effort to get the best fit 
to the data. The coupling constant was determined for 
each mixing angle from the given total width of the a\ 
resonance. Besides the intermediate states with the a\ 
resonance, we considered also those with only pions and 
p's as given in various theoretical schemes [!, H, H, [l(| • 
They influence the calculated cross section mainly in the 
rho mass region. However, the quality of the fit through- 
out the energy region covered by BaBar experiment [131 ] 
has improved only slightly. 

As far as the e + e~ annihilation cross section data are 
concerned, the situation in the 7r + 7r _ 7r°7r° sector is worse 
than in the 7r + 7r~7r + 7r _ one. The data coming from var- 
ious experimental groups did not agree with one another 
very well; see, e.g., Fig. 10 in 14]. It is good news 
that the latest published data by the SND Collaboration 
at BINP in Novosibirsk agree well with the newest 




FIG. 1: Two Feynman diagrams of a pure-ai model of 
e + e~ — ■» 7r + 7r~7r°7r°. Two others can be obtained by exchang- 
ing 7T°'S. 

BaBar [3 [H| and CMD-2 [3 [l7| data. Unfortunately, 
those data are still preliminary and publicly unavailable 
.18, 19]. We can therefore use only the SND data, which 
cover the energy interval from 980 to 1380 MeV. The sta- 
tistical errors are combined with the 8% systematic error 
[Til ] in quadrature. 

We will first compare the cross section data to a simple 
pure-ai model, which is characterized by four Feynman 
diagrams depicted in Fig. [1] The model is an obvious 
modification of the model used in [7j , which is defined by 
a two-component aipir interaction Lagrangian 

C ai(m = ^-(C lC ose + C2sme), (1) 

where 

A = A" • (V^ x d v <j>) , (2) 

C 2 = • (d»A v x 0) , (3) 

and V M „ = d^Yu — <3„V M . The isovector composed of the 
/3-meson field operators is denoted by V M ; similar objects 
for 7T and a\ are <j> and A**, respectively. The sine of the 
mixing angle sin (9 = 0.4603(28) was determined in Q 
by fitting the e + e~ — * tt + tt^tt + tt~ cross section data 
from the BaBar Collaboration fl3[ supplemented with 
the experimental value of the D/S ratio in the a\ — > 
pn decay [2fj| . The value of the coupling constant g aiP Tr 
follows from sin 9 and the total width of the a\ meson, 
which was chosen at 600 MeV. 
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FIG. 2: The generic Feynman diagram describing the ui con- 
tribution to e + e~ — > 7r + 7r~7r°7r . The other five diagrams are 
obtained by obvious modifications. 



FIG. 3: One of the two Feynman diagrams with the contact 
u>3tt term. Another is obtained by exchanging 7r°'s. 



Also defining our model is the form factor generated by 
the p(770), p' = p(1450), and p" = p(1700) resonances, 



F(s) = F p (s) + 6F p ,(s)+eF p „(s). 



(4) 



As far as the individual contributions on the right-hand 
side are concerned, we refer the reader to formulas in 
Q. Here we only note that the complex parameters e 
and S as well as the masses and widths of the p' and p" 
resonances hidden in F„i and F p n were determined by 
fitting the four-charged-pion BaBar data |13j . 

The last ingredient of our model is connected with the 
structure of the strongly interacting particles. Each ver- 
tex is usually modified by a strong form factor to soften 
the interaction. In 0], we used a simplified approach. 
We merged all form factors to one, effective, strong form 
factor of the Kokoski-Isgur 21] type, which multiplies 
the total annihilation amplitude 



Fki(s) = exp <^ - 



s - s 
48/3 2 



(5) 



where s = 16m 2 and (3 = 0.3695(98) GeV follows from 
the fit to the BaBar data [llj]. The same form factor is 
also used here. 

When calculating the e + e~ — * 7r + 7r~7r 7r° excitation 
curve in the pure-ai model, we keep all the parameters 
at values determined in Q. The result \ 2 = 2076 for 
35 data points is disastrous. A poor result for the pure- 
ai model clearly signifies that an additional contribu- 
tion to the amplitude of the electron-positron annihila- 
tion into two charged and two neutral pions is needed. 
The intermediate states with the to meson, considered al- 
ready by Renard in 1969 [l[ and later by other authors 
[1,0, SH, [H|, are an obvious choice. 

We will pursue two different ways of including the in- 
termediate states with the u> meson. First, we adopt the 
approach of Eidelman, Silagadze, and Kuraev (ESK) Q, 
who used the anomalous part of the chiral Lagrangian 
[22l l23l |24| to describe the ivpir and u)3ir vertices. The 
Feynman diagrams are shown in Figs. [5] and [3J The 
second approach (PL) 25] is more phenomenological. It 
does not consider, like [3j,[ll|, the w37r contact term. The 
Lagrangian 



G 



(0 ■ d p V a ) 



(6) 



has the same form as in 0, 0, [Tl| . But now, the coupling 

37r decay 



width taking into account the value of the pirn coupling 
constant g p as it follows from the p — > -kit decay width, 
g 2 p = 35.77 ± 0.24. We get G 2 Ufm = (216.2 ± 3.0) GeV" 2 . 
The value of G^ p7T is higher than the corresponding quan- 
tity in [9j by only about 2.6%, so the main difference 
between the two approaches lies in the diagrams with 
the uj3tt contact terms. To check the soundness of our 
approach, we performed two tests. 

First, we calculated the width of the radiative decay 
lo — > 7r°7 assuming the strength of the p°j coupling, as 
it follows from the normalization of the pion form factor, 



em 

£"f P ° = 

9 P 



( 7 ) 



constant G wp7r is determined from the ui 



The calculated branching fraction of (9.48 ± 0.28)% dif- 
fers a little from the current experimental value B{lo — > 
^ 7 ) = (8.90i° ;^)% M- 

Second, we determined the strength of the lvj cou- 
pling from the w — > e + e~ decay width and used it in 
calculating the rate of the 7r° — > 77 decay. The re- 
sult, expressed in terms of the tt° mean lifetime, t = 
(7.7 ± 0.4) x 10~ 17 s agrees well with the experimental 
value of (8.4 ± 0.6) x 10~ 17 s. 

What remains unsettled is a possible transfer- 
momentum-squared (t) dependence of the pj coupling. 
In fact, the experimental width of the p° — > e + e~ de- 
cay, where t — m 2 pl requires about 20% stronger coupling 
than that indicated in ([7]). The latter gives good results 
for the t = processes lo — > ir°j and 7r° — * 77. Our 
derivation Q of the e + e _ — > 47r cross section formula 
used the standard 7/5 coupling ([7]) and assumed that all 
the t dependence is absorbed in the form factor. We use 
the same approach here. 

Now, we can also vary the form-factor parameters 5 
and e, as the structure of the intermediate states is dif- 
ferent from the pure-ai model. We will distinguish them 
from the 7r + 7r~7r + 7r~ case by primes. There is no free pa- 
rameter connected with the lot: intermediate states. The 
results are shown in Table U 

In an effort to further improve the agreement of our 
model with data, we include the intermediate states with 
the isoscalar axial-vector meson /ii(1170). The corre- 
sponding Feynman diagrams can be obtained from those 
in Fig. [2] by replacing w with h\ . The interaction La- 
grangian is again assumed in a two-component form 
similar to (Q]) but respecting the isoscalar character of 
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TABLE I: Fitting the ax+uo model to the e + e~ 


-> tt+tt-ttV 


cross section data (35 data 


points). 




Approach to to 


"c 1 c T/" rnl 

EbK [^| 


PL [25] 


v 2 /NDF 


2.19 


0.82 


CL (%) 


0.01 


7 A Q 


Re 5' 


-0.52(25) 


-0.36(27) 


Im S' 


-1.27(25) 


-1.05(32) 


Re e' 


-0.79(34) 


-0.71(37) 


Im e 


0.953(97) 


0.628(79) 


TABLE II: Fitting the ai+u+hi model to 


the e"^e~ — ► 


7r + 7r _ 7r°7r° cross section data (35 data points) 




Approach to ui 


ESK [9] 


PL [25] 


X 2 /NDF 


1.18 


0.80 


CL (%) 


22.4 


77.8 


sin rj 


0.3434(36) 


0.3433(46) 


Re 5' 


0.092(62) 


0.102(74) 


Im 5' 


0.028(22) 


0.035(23) 


Re e' 


0.022(71) 


0.028(86) 


Im e 


-0.030(58) 


-0.049(65) 



/ii(1170), 

^fe lP 7r = g ^- ?r (£a cos 77 + C b sin??) , (8) 

where 

C a = ^-(V^.9», (9) 
C b = -W-(V^-<^). (10) 

In the following, the sine of the mixing angle rj will be 
varied to achieve the best possible description of the cross 
section data. For each sin 77 the coupling constant will be 
determined from the condition that the total width of the 
/ii(1170), calculated as T(h\ — > pit), should be equal to 
360 MeV. While fitting the e+e~ -> 7r + 7r^7r 7r data the 
mixing angle of the aipn Lagrangian ([1]) is kept fixed at 
the value determined in the 7r + 7r~7r + 7r~ case, as it rep- 
resents a universal process-independent parameter. The 
results of the fit are shown in Table HT1 for both approaches 
to the intermediate states with lu. It is clear that the in- 
clusion of the hiir intermediate states greatly improves 
the confidence level of the model with ui described by the 
ESK scheme. The confidence level of the model utilizing 
the PL scheme for lu also rises, but because it was already 
high in the cli+lu model, the inclusion of the intermedi- 
ate states with hi is not necessary. The excitation curves 
are compared to data in Fig. |U 

Following the idea that a simultaneous fit to more pro- 
cesses may lead to a more precise value of the mixing 
angle of the aipir Lagrangian, we are now going to per- 
form a joint fit to the e + e~ — > 7r + 7r~7r + 7r _ and e + e~ — > 




W (GeV) 



FIG. 4: Comparison of the ai+u)+hi model with the e + e — > 
7r + 7r - 7r°7r° data for two approaches to lu. 

7r + 7r~7r°7r° cross section data. We merge the pure-ai 
model of our previous paper Q with the ai+uj+hi model 
described here. For the ujtt intermediate states we will 
use the PL version, which describes the data better than 
ESK. Simultaneous handling of both four-pion channels 
enables us to use a more correct description of the p(770) 
part of the electromagnetic form factor Q , namely, 

M„ 2 (0) 

FJs) = pV ' . , (11) 

P M*{s)-s-im p T p {sy 

where M (s) is the running mass of the p meson calcu- 
lated in [27| • The total decay width of the p° meson, 

r p (s) = T p (s) + r 7r+7r - 7r+Tr - (s) + r^+^-^o (s) (12) 

now includes not only the contribution of several two- 
and three-body decay channels T p (s) given in [27J, but 
also the contributions from both four-pion decay modes. 
Similarly to we also consider the masses and widths 
of p 1 and p" as free parameters. They have the same 
values in both channels, as well as the aipir Lagrangian 
mixing parameter sin 9 and the parameter /3 of the strong 
form factor ([5|). Besides those six common parameters, 
there are two sets of parameters specific for each of the 
two four-pion annihilation channels. The charged-pion- 
channel set contains four real parameters entering the 
electromagnetic form factor The four form- factor 

parameters in the mixed-pion channel are distinguished 
from those in the 7r + 7r~7r + 7r~ channel by primes. The 
fifth parameter in the mixed-pion channel is the h\pix 
Lagrangian mixing parameter sin 77. All together, this 
makes 15 real free parameters. 
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TABLE III: Joint fit of the pure-ai model of e + e~ — > 
Tv + ir~ir + iT~ and the ai+tj(PL)+hi model of e + e~ — > 
7r + 7r - 7r°7r° to a combined set of the cross section data and 
the D/S ratio (180 data points). 



Quantity 


Value 


Quantity 


Value 


X 2 /NDF 


1.06 


Re 8' 


0.088(73) 


77V ( GeV ) 


1.383(16) 


Im 8' 


0.024(22) 


r> (GeV) 


0.551(21) 


Re e' 


0.025(81) 


m p „ (GeV) 


1.883(18) 


Im e' 


-0.011(74) 


F p „ (GeV) 


0.237(35) 


Re <5 


0.171(14) 


sin 6 


0.4662(52) 


Im <5 


0.045(29) 


(3 (GeV) 


0.3617(70) 


Re e 


0.0002(14) 


sin r\ 


0.336(11) 


Im e 


-0.0050(12) 



In Tabic IIIII we present the optimized values of all 
free parameters. The corresponding confidence level is 
28.4 %. If we compare them with those obtained in in- 
dividual models (Table HT1 here and Table VI in Q), we 
can say that they are in good agreement. Only the er- 
ror domains of m p i, m p ", and d do not overlap, but the 
disagreement is very small. The excitation curves of in- 



dividual reactions do not visually differ from those ob- 
tained when fitting the two models separately and are 
not shown. 

Unfortunately, our goal to narrow the interval of the 
aip-K Lagrangian mixing parameter and thus make the 
calculations of the dilepton and photon production from 
hadron gas more reliable (see the analysis in [28[) has not 
been reached. The uncertainty of sin 9 is larger than that 
obtained in [7j. We hope that the situation will improve 
when more precise cross section data in the mixed-pion 
channel arc available. Identifying the essential contribu- 
tions to the electron-positron annihilation into the four- 
pion final states is important for the reliable assessment 
of the dilepton production by the four-pion annihilation 
in hadron gas [291 ] . A new result of our work is the mixing 
angle of the hipir Lagrangian sin 77 ss 0.34. It may help in 
investigating the role of the ft.i(1170) resonance in ther- 
mal production of dileptons and photons from hadron 
gas, which has been ignored so far. 
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